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Abstract. Collisionless shocks occur in various fields of physics. In the context of 
space and astrophysics they have been investigated for many decades. However, a 
thorough understanding of shock formation and particle acceleration is still missing. 
Collisionless shocks can be distinguished into electromagnetic and electrostatic shocks. 
Electromagnetic shocks are of importance mainly in astrophysical environments and 
they are mediated by the Weibel or filamentation instability. In such shocks, 
charged particles gain energy by diffusive shock acceleration. Electrostatic shocks are 
characterized by a strong electrostatic field, which leads to electron trapping. Ions are 
accelerated by reflection from the electrostatic potential. Shock formation and particle 
acceleration will be discussed in theory and simulations. 


1. Introduction 

Collisionless shocks show interesting features for particle acceleration, which is why they 
are of importance in many helds of physics ID-IZl- They are generated by the interaction 
of charged particles with the surrounding self-generated helds - in contrast to collision- 
dominated shocks, where two-particle interactions determine the physical behaviour. 
In a simple setup of two interpenetrating plasma slabs, either a strong electrostatic 
held [8] is generated due to the two-stream instability [9], or an electromagnetic held 
due to Weibel-like instabilities mm- The feedback of such helds mediates the shock 
formation process. 


2 


The Weibel or filamentation instability generate a turbulent field in which the 
charged particles are defiected. Parallel momentum is transferred into perpendicular 
momentum, and thus the particles are accumulated and the plasma slabs are compressed. 
In relativistic initially unmagnetised plasmas, the density compression is three times the 
initial density HZ], The time scales of the shock formation are determined by the time 
scales of the electromagnetic instability [13] . The details are given in Sec. [3l 

Electrostatic (ES) shocks form in plasmas with different components of electrons 
and (heavy) ions. A high mass and temperature difference is favourable for the 
generation of these shocks. An electrostatic potential is built-up which traps electrons 
in the downstream region and the electron distribution function is widened in parallel 
direction. The perpendicular components are less affected. This temperature anisotropy 
gives rise to the electromagnetic Weibel instability. The time scales of the formation of 
such electromagnetic modes are calculated in Sec.lHand compared against the time scales 
of shock formation in order to determine their relevance during the shock formation 
process. We perform the analysis for plasmas with relativistic temperatures and fluid 
velocities. Laser-generated plasmas in the laboratory are now entering the relativistic 
regime, which is why it is important to explore the physics also in this parameter range. 

2. Initial setup for shock formation 

We study shock formation in a simple symmetric setup (see Fig. [T]). Let us consider 
two charge-neutral counterstreaming beams of electrons and positrons or ions with 
bulk velocities ±uo, temperatures Tg, Tj and thermal parameter /i = mc^/ksT. The 
overlapping region turns unstable due to collisionless plasma instabilities (Phase 1), 
which will mediate two shocks propagating into the upstream regions (Phase 2). 


Phase 1 


Unstable 

region 




Figure 1. Shock formation from the interaction of two counterstreaming beams. 
In Phase 1 the overlapping region becomes unstable. In Phase 2 two shocks are 
propagating into the upstream regions. 
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The collisionless shock is of electromagnetic or electrostatic nature, depending on 
the initial parameters for the velocities and temperatures. We will discuss these details 
in the following sections. 


3. Shock formation in electromagnetic shocks 


Electromagnetic shocks develop if electromagnetic instabilities dominate in the 
overlapping region. In such a symmetric setup, the electromagnetic hlamentation 
(Weibel) instability is the fastest mode for cold beams Te ~ T* 0 and relativistic 
velocities v^jc ~ 1. The growth rate of the cold hlamentation instability is given by 
= ^^ 2 / 70/90 ‘^pa with /3o = Vq/c, Lorentz factor 70 = (1 —and plasma frequency 
of species a being ojpa = Annoo^/rria. The saturation time of this instability 


Bj 


2Sa 


( 1 ) 


is a function of the initial and hnal magnetic held strengths. The initial magnetic held 
strength Bi is obtained from the evaluation of the spectra of spontaneous magnetic 
huctuations, while the hnal held 5/ — Sir'yQnomc^ results from a trapping condition 
in the magnetic held structure mi. The saturation time in electron-positron shocks can 
then be expressed as 

3 \ 


^s,e^pe 





( 2 ) 




For electron-positron pair shocks, the shock formation time is simply twice the saturation 
time of the instability in two dimensions, Tf^e = 2Ts^e, and a factor 3 in three dimensions, 
rf,e = US]. 

For electron-ion shocks the theory has to be extended by an extra term for the 
merging time of the hlaments. At the saturation time of the hlamentation modes in 
electron-positron plasmas, the transverse size of the hlaments is already large enough 
in order to dehect particles strongly enough for efficient accumulation of particles. On 
the other hand, in electron-ion plasmas, the hlaments are still on the electron scale. An 
additional merging time 

93/2 

. ( 3 ) 
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is required in order to bring the hlaments to the required size to signihcantly dehect the 
ions [ 12 ]. The total shock formation time in electron-ion shocks is thus given by 

Tf,i ojpi = {Ts,i + Tm) ojpi = 4.43 d 7 ^^ ln(mi/ rrie) (4) 


with Ts^i the saturation time of the cold ion hlamentation instability and d the number of 
dimensions. The predicted scaling is broadly consistent with particle-in-cell simulation 
results. 
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4. Shock formation in electrostatic shocks 


Electrostatic shocks require a mass and temperature difference between ions and 
electrons. The electrostatic two-stream instability is dominant in the overlapping region, 
which is fast for low streaming velocities. A strong electrostatic potential is generated in 
the downstream region which can trap electrons (see Fig. [2]), which can be approximated 
as a flat-top distribution. In the non-relativistic case, we follow Ref. [16] for a steady- 
state solution of the shock and describe the electron distribution by a population of 
free streaming and a population of trapped electrons. If we introduce normalisations 
for the velocity, /3 = n/c, and the electrostatic potential, = e^/mgC^, electrons are 
free streaming if the condition \l3x\ > is fulfilled and trapped otherwise. The non- 
relativistic distribution of the electrons is then given by [16] 


fe 


Co 


2n) 


exp{-/i(\//?2 _ 2^ + Pof/2} 

exp{-/r/lo/2} 

exp{-/i(y//32-2(/?- (3o)^/2} 


f^x ^ \/ 2 (^ 

1/5x1 < ^/2^ (5) 
/5x ^ \/2^ 


with the normalisation constant Co = 


e^'^ericy/JIip + 2y^/i^7^e 


-1 


Fig. El 


demonstrates the ES shock configuration with an oscillatory potential in the shock 
downstream region. 



Figure 2. An electrostatic shock forms in the interaction region of two counter 
propagating plasma slabs. The electrostatic potential y increases monotonously from 0 
in the upstream (a), a monotonous increase in the transition region (b) to its maximum 
‘fimax in the downstream (c). 


4-1. Relativistic generalisation of the distribution function 

For the relativistic generalisation of the electron distribution eq. (jS]), we start with a 
Maxwell-Jiittner distribution in the mean rest frame (subscript R) 

f^^ = Cr exp (-/iij 7 it), ( 6 ) 

with the Lorentz factor always defined as 7 = (1 — and perform a Lorentz 

transformation with 70 into the moving (= laboratory) frame (subscript L) in ±x 
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directions. Thus, with 7 ^^ = 7o7l(1 ± Pof^L,x) we obtain the relativistic distribution 
in the laboratory frame 

= Cl exp (-/iL 7 L (1 ± ^o/3l,x)) (7) 

with fiL = hJ? 7 o- Cr and Cl are normalisation constants in the mean rest frame and in 
the laboratory frame, respectively. 

We now introduce the dependency on the electrostatic potential (p. From energy 
conservation, we obtain a balance between the kinetic energy in the upstream, where 
the bulk is moving with 7 /,, and the total energy in the shock transition region, where 
kinetic energy is transformed into potential energy. The bulk Lorentz factor is reduced 
to 7 ^. This is expressed by 


meC^i'yL - 1 ) = meC^i'JL - 1 ) - 60 


( 8 ) 


or simply Jl = Jl ~ 9 ^- Electrostatic shocks are mainly one-dimensional, which is why 
we assume that the longitudinal and transverse processes can be separated. Introducing 


1 . 


the dimensionless momentum u = 07 , we approximate ul,x ^ ^ '^'1 x ~ 

We can then write the relativistic generalisation of the electron distribution in eq. (1^ . 
if we replace ■jl and (3l,x in eq. ([7]) by their ip dependent terms and with 7 := 7 ^ and 


Ux := a, we obtain 


exp < —/i 


7o (7 - - 1 + “0 


07 i+ 


ut 




fre{u) = CrO exp {-/i [7o7± - 1]} 


exp < —p 


7o 


(7 - 9 ?) - 1 - uo\j+ ipY - I 


with the dehnitions 7 = 7 (u) = a/I + and 7 _l = + u\. The dehnition of the 

trapping velocity Uc is derived from the balance of electrostatic and kinetic energy in 
eq. dH]) for 72 , = 1 which gives 7 = 7 c := 1 -|- y? and Uc = \/'il — 1. The non-relativistic 
approximation of this trapping velocity agrees with the non-relativistic dehnition in eq. 
dS]) for Me — )■ 0c = a/2^ for 0c <C 1. 

The normalisation constant Cro is obtained from 


d^u fr. 


lu, 


= 1 


( 10 ) 


with 


CrO — 


27re^ 


Ux 


m, 


Ux 


2mc( 1 + 7oh)e"'^'^° + ^ 


'9c 


T 7 ^ 707 ) ^[-;j(707±M0^/(7-y)^ 


Eqs. dS]) and dS]) describe the electron distribution in the quasi-steady ES shock. 
In hg. [3] we show the change of the distribution function from the initial Maxwell 
distribution for 99 = 0 to a hat-top distribution for 99 > 0 for uq = 0 o 7 o = 0.01 
and p = 50, corresponding to ksTe = 10 keV. The relativistic expression eq. (jH]) is 
also compared against the non-relativistic expression eq. (El)- In the far upstream. 


< -Uc 

I < Uc (9) 

> Uc 


-‘>1 ( 11 ) 
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the distribution is Maxwellian, see fig. In the shock transition region, where 

0 < < ipmax^ the distribution f{ux) broadens and becomes fiat-top. (pmax denotes 

the saturation value of the electrostatic potential in the far downstream which can be 
obtained by solving numerically the differential equation | = 0 with the 

Sagdeev potential [IZl HE]- Fig. [3b shows the distribution function at the ion 

reflection condition, where the kinetic energy of the ions equals the electrostatic energy. 
In normalised quantities, this is given hy ip = (prefi •= (7o ~ In the far 

downstream, for ip = (pmax the non-relativistic description breaks since (pmax > 1- 




Figure 3. Electron distributions for ip = 0 (a), (prefi = 0.09 (b) and Pmax = 4.58 (c) 
for uq = 0.01 and p = 50 obtained from the relativistic expression eq. (|9]) given in 
blue and the non-relativistic approximation eq. ([5|) in dashed black. 


In this configuration, the electrons are usually treated kinetically, while the higher 
inertia ions are described with a fluid model. In the following section, we summarise 
the initial conditions for which an ES shock forms. 


4-2. Conditions for ES shock formation 


Here, we summarise the results on ES shock formation conditions. In Refs. mm 
a condition for the maximum Mach number was found as 1 < Mmax ^ 3.1, with 
the Mach number M defined as the ratio of the upstream velocity in the shock rest 
frame to the ion sound speed, M = Vq/cs = (uq + Vsh)/cs. The shock velocity in 
the laboratory frame can be approximated from the shock jump conditions as Vsh/c = 
■s/Sad — l\/(7o “ l)/(7o + 1) with the ideal gas adiabatic constant yad [I2]. The ion 
sound speed is given by Cg = ^yksTe/mi. This can be generalised for relativistic plasmas 
with the relativistic Mach number JCl = Uq/us, where Wq = /9o7o is the dimensionless 
momentum in the shock rest frame and Ug = /9s7s with (4g = ^^ksTe/rUi c‘^ mm- This 
imposes a condition for the upstream fluid velocity and electron temperature 



1 < < 3.1. 


( 12 ) 


5. Calculation of unstable modes 


During the early stage of ES shock formation the broadening of the electron distribution, 
which we observed in the previous section, occurs mainly in the longitudinal direction. 
The transverse directions stay almost unaffected. This has also been observed by 
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particle-in-cell simulations [22]. The generated temperature anisotropy in the electron 
distribution gives rise to electromagnetic modes. 

We develop now a model to describe the growth rate of the electromagnetic modes 
in such a setup. Later, we will compare them to the shock formation time scales and 
the growth rate of the cold ion-ion instability in order to determine their relevance. 


5.1. Dispersion relation of EM waves in relativistie plasmas 

We start from the electron distribution in eq. ([9|) in order to evaluate the dispersion 
relation of electromagnetic waves in plasma 

fcV-a;2 

with 

poo 

u = 


and 


1 / = 


- + r) = 0, 

(13) 

,3 Ux df 
d-^u 

7 

(14) 

<eu 

(15) 


We consider only fluctuations k = ke^ perpendicular to the fluid velocity uq = MoGx in 
order to simplify the geometry. An evaluation of the integrals leads to 

^ f y'y _ XgwV(T'-¥>)^-i 

± “'7c 

(3o 7 - 9 ? 


U = 


X 


7o7r (0, /i7o7) ± 


„-^707 


y(7 - - 1 


(16) 


and 


V = Cro4^7o 


du^u^ 


dUn, 


3-M(707x-1) 


dUn 




,7x 


7(7V + ul) 


E- 




ut 


72(72^2 ^.^2) 


3-M(707±«0^(\/l+Mi-¥’)^-l) 


(17) 


where the last integrations have to be done numerically for the general case. 

We use eqs. flTbl) and ffT7|) to solve the dispersion relation flT3l) numerically to obtain 
the growth rate <j{k) := Q{u{k)). In Figs. and b we plot the maximum growth rate 
o'max '■= raax{a{k)) for different values of p and Uq. The maximum growth rate was 
evaluated at the ion reflection condition (frefi = (7o~l)n^i/''ne. The role of the potential 
will be discussed in sec. 15.31 

o'max shows a maximum as a function of the velocity. The higher the electron 
temperature, i.e. the lower /i (x the higher the velocity at which the maximum 
of O'max appears (£g. |4^). On the other hand, Omax increases with /i. For an easier 
interpretation, we plotted the dependence of Omax against the electron temperature in 
Fig. IDd. The higher the initial electron temperature, the lower is the growth rate Omax- 
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Figure 4. Maximum growth rate amax vs. uq (a) and fcsTe (b) obtained from 
dispersion relation eq. m with relativistic expressions eqs. (USD and (II3- 


5.2. Dispersion relation of EM waves in non-relativistic plasmas 


In order to further discuss the properties of this model and to check the consistency with 
previous models, dispersion relation flTdll is approximated for non-relativistic velocities, 
i.e. /So ^ 1 and 1. We obtain Un ~ —1 and 


Vn-Co 






e'^'^erfc-y/p^ + 




■M/3o/2 


.(18) 


For zero beam velocity and </? = 0, this is exactly the well-known solution of the Weibel 
instability m- For small frequencies u <C kc, the plasma dispersion function Z can be 
further approximated and we get 




1 






eT^c^/JIip + 2 




7^So/2 


.(19) 


The dispersion relation then reads 


7 2 2 2 I 

k C — OJ + OJ 


2 

pe 


1 - V{^) 



0 


( 20 ) 


with V(cp) = Co< e^'^eiic^/Jup + 


+ 


; /ZZe-M/3oV2 

1 y TT 


The solution of the 


dispersion relation can now be derived analytically, which is given by 


a{k) = 2s{uj{k)) 



kc 


/ivr 




( 21 ) 


with k^c^ = ujp^{V{(p) — l)/3 the location of the maximum and the maximum growth 
rate 


1)) . (22) 

For comparison, fig.[5]shows the relativistically correct solution of the dispersion relation, 
obtained from eqs. ffTOj) and flTTIi . and the non-relativistic approximations from eqs. fl 20 |) 
and (El]). The parameters used are (3o = 0.01, p = 100 and (p = 0.05 (70 —1) mi/me. The 
maximum growth rate matches with the approximation in eq. fl 22 |) . (Jmax = 2.0 x lO^^Wpe 
and the location of the maximum at fco = 0.25a;pe/c. 
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Figure 5. Numerical solution of the dispersion relation for the relativistically correct 
eq. (1181) in blue, the non-relativistic approximations from eq. ( 1201 ) in black and from 
eq. (1^ in red for /?o = 0.01, ^ = 100 and ip = 0.05 (70 — 1) mi/me. 


5.3. The role of the electrostatic potential 

The electrostatic potential (p is quickly built up during the electrostatic shock formation. 
Nevertheless, the electromagnetic growth rate a{k) depends on (p and thus, the 
dependence should be investigated. We calculate the growth rate numerically for values 
0 < </? < (prefi, which are plotted in £g. [6l 



Figure 6 . Maximum growth rate for uo = 0.01 and /r = 100 (yellow), 10 (orange), 
1 (red), 0.1 (brown) and uq = 0.1 and p = 1 (light blue), 0.1 (dark blue) with 

^ ^ [ 9 , Prefl \■ 


We hnd a power-law dependence of the growth rate on the electrostatic potential, 
which gives 

^max T (^ 2 ) 

with a > 0. This ip dependence can be interpreted as the changing growth rate across 
the steady-state shock. If we assume, as a simple approach, that the potential grows 
linear in time, we can connect the growth rate ccmax with time, which will be further 
discussed in the next section. 
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6. Comparison of different time scales 

In order to determine the relevance of the electromagnetic electron Weibel modes during 
the ES shock formation process, we compare the growth rate given by eq. fll3l) . which 
we label from now on aEM,ee, with other processes. 


6.1. ES Shock formation time 

The time scales of ES shock formation are determined by the growth rate of the ion-ion 
electrostatic instability which is given by [9] 

0 'ES,ii = -372 

27o 

Estimating the shock formation time with 5 times the inverse maximum growth rate, 
we obtain tgf = mi/meUpf, which was conhrmed in simulations (see [23]). 


6.2. Growth rate of the electromagnetic cold ion-ion instability 


Another competing process in ES shocks is the cold ion-ion hlamentation instability, 
which has a growth rate [ 23 ] 



(25) 


O'EM,a — (do 


The EM mode of the ion-ion instability grows faster than the ES mode for 
fluid velocities uq/c > 1/3. In the parameter range vq/c < 0.1, which we looked on in 
this paper, the ion-ion EM modes can be neglected. 


6.3. Dominant regimes 

A comparison of the ES shock formation time scales with the EM modes makes it 
possible to determine parameter regimes, for which a shock stays electrostatic. For 
this, we compare the growth rate from eq. (fT3|) with eq. (|24|) for different electron 
temperatures, expressed by /c^Te/meC^, and fluid velocities, given by Uq = /So7o- For 
simplicity, we choose an electrostatic potential (p = p>refi- 

In the dark green region in Eg. [7] the growth rate of the EM electron instability 
OEM,ee E Smaller than the ES ion-ion growth rate OES,ii- We label this domain the 
purely ES domain. In contrast, in the light green region, an ES shock will turn EM 
since OEs,ii < OEM,ee- In the white region (EM), no electrostatic shock develops because 
condition ffT2|) is not fulhlled. This regime is electromagnetically dominated from the 
beginning, with shock formation according to section [S] 

A more detailed analysis of the scenario in Eg. [7] has been done in Eg. [H] for difierent 
values of the electrostatic potential (p. For difierent fluid velocities Uq and temperature 
parameters p = meC^/ksTe, the maximum growth rate amax is calculated as in fig. [6l 
The characteristic time is then calculated as tchar ~ ^(omax and plotted against the 
time of shock formation, where we assumed that the potential grows linearly as (p (x t. 
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Figure 7. Determination of dominant regimes as a function of dimensionless 
momentum u = /3o7o and electron temperature in keV, cf. Ref. [23]. 



Figure 8. Characteristic time tchar ~ ^Ic^max against time during shock formation 
t € [0, 107 q with same colour coding as in fig. |6| The grey region shows 

ichar < C where EM inodes are important during ES shock formation. 


The grey area in fig. [HI shows the region, where EM modes develop faster than the 
ES shock, which is the case e.g. for mq = 0.01 for /x > 10 (< 51 keV) or for uq = 0.1 
for /i > 0.1 (< 5 MeV). Fig. IHlalso shows that a detailed analysis of the growth rate is 
not necessary, since the characteristic time tchar quickly saturates. The rough estimate 
in fig. |6]is thus sufficient for a qualitative determination of the regimes. 

7. Summary and conclusions 

Electromagnetic and electrostatic shocks can both develop from the same simple 
symmetric setup of counterstreaming beams. The choice of the initial plasma parameters 
determines the final shock character. 

Electromagnetic shocks develop in plasmas with large beam velocities and low 











12 


temperatures where the electromagnetic filamentation instability is the fastest mode. In 
electron-positron pair plasmas, the shock formation time is simply twice the saturation 
time of the instability in 2D. In electron-ion plasmas an extra merging time is necessary 
in order to obtain the right scale of the filaments to efficiently scatter the ions. Thus, 
the shock formation in electron-ion plasmas is delayed by almost a factor 3 compared 
with the pair shock, Tf^iOJpi ~ 3r/_e<^pe- 

Electrostatic shocks form in electron-ion plasmas with small beam velocities and 
large electron temperatures. Here, the electrostatic two-stream instability dominates. 
An electrostatic potential is generated in which electrons are trapped. The subsequent 
deformation of the electron distribution gives rise to electromagnetic Weibel modes 

which can destroy the electrostatic shock features. 
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